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The decline in extent of wild pollinators in recent years has been partly associated with 
changing farm practices and in particular with increasing pesticide use. In this paper we 
combine ecological modelling with economic analysis of a single farm output under the as¬ 
sumption that both pollination and pest control are essential inputs. We show that the drive 
to increase farm output can lead to a local decline in the wild bee population. Commercial 
bees are often considered an alternative to wild pollinators, but we show that their intro¬ 
duction can lead to further decline and finally local extinction of wild bees. The transitions 
between different outcomes are characterised by threshold behaviour and are potentially 
difficult to predict and detect in advance. Small changes in economic parameters (input 
prices) and ecological parameters (wild bees carrying capacity and effect of pesticides on 
bees) can move the economic-ecological system beyond the extinction threshold. We also 
show that increasing the pesticide price or decreasing the commercial bee price might lead 
to re-establishment of wild bees following their local extinction. Thus, we demonstrate the 
importance of combining ecological modelling with economics to study the provision of 
ecosystem services and to inform sustainable management of ecosystem service providers. 
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I. INTRODUCTION 


Globally, 35% of food crops are at least partly dependent on insect pollination ([Klein et ^ 


20071. Ensuring sufficient pollination of these erops will be ehallenging in the future; the fraetion 


of agriculture made up by inseet-pollinated crops is increasing ( jAizen & Harder 2009), while wild 
pollinator populations are threatened by both habitat loss ( jWinfree et ah 20091 and agricultural 
intensification, which are thought to be the main eauses of reported deelines in diversity in the EU 
and in the USA ( jBiesmeijer et al.\ 2006; Cameron et al. \ 201 Ij ). 

Eor some crops, honeybees are used to supplement or substitute wild pollina¬ 
tors, along with other commereial pollinators sueh as laboratory bred bumblebees 
( Velthuis & van Doorn[ 20061. While eommercial pollinators are often assumed to be adequate 
substitutes for wild pollinators (though see ( [Brittain et al.\ |2013[ [Hoehn eiW] |2008[ )), the use of 
commereial pollinators is itself not without risk. Honeybees have suffered losses in reeent years 


due to the abandonment of hives (Colony Collapse Disorder) and the Varroa mite (Cox-Eoster 


et al. 2007). Relying on eommercial pollinators such as honeybees puts farmers at risk from 


shocks of this kind, with consequent implieations for farm profits over time. 

Given the risks around the supply of pollination serviees from commereial bees, maintaining 
viable wild pollinator populations is likely to be crueial to sustaining the produetion of insect- 


pollinated erops into the future ( [Winfree et al. \ [2007 p . The potential eosts of the loss of loeal 
pollination serviees is illustrated by the need for pollination using farm workers in Sichuan, China, 


following the loss of loeal insect pollinator populations (Partap et al.. 20011. Whilst this was a 
viable option when wages were eheap, a 10-fold rise in wages over the last 10 years has led to 
the abandonment of apple produetion ( Partap & Tang[ 2012). One of the factors implicated in this 
loeal extinction and in declines elsewhere, is the use of pestieides, or specifieally, inseeticides. 
There is growing evidence of negative effeets of realistic levels of eommonly used insecticides 
on population determining traits sueh as reproductive rates, foraging rates and navigation in bees 
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(Goulson 2013[ Henry et al. 2012 Mommaerts et al. 2010[ Whitehorn et al. 20121. Awareness 
of this evidenee has led to the temporary banning of a very widely used group of inseetieides - 
neonicotinoids - within the European Union. 

Farmers who grow inseet-pollination dependent crops face a trade-off in their use of pesticides, 
since whilst this reduces crop damages, it also has potential negative effects on the local supply 
of pollination services from both wild and commercial pollinators. This is the issue that we study 
in this paper. To investigate the links between commercial pollinator use, pesticide use and wild 
pollinator populations we present an ecological-economic model which links crop yields to polli¬ 
nator numbers and pesticide use, and study how the optimal strategy depends on the level of farm 
output and other parameters. 

There are a number of key results in this paper: Firstly, by introducing commercial bees the 
farmer can achieve higher target output values than by relying on wild pollinators alone for the 
same level of cost. Thus, it is rational to substitute commercial for wild bees and to increase 
the use of insecticides. Secondly, the resulting increase in pesticide use buffered by usage of 
commercial bees may lead to a severe reduction in wild pollinator population, even without direct 
competition between different bee populations. Thirdly, under certain economic conditions, the 
system becomes unstable and small changes in prices or environmental conditions can result in 
the local extinction of wild pollinators. As in this paper we are modelling a single farm only, all 
references to extinction mean that locally the wild bee population will tend to zero. 

Finally, we discuss a number of options that are available to restore a wild bee population even 
if it is currently locally extinct. These include an increase in the pesticide price (or application 
of a pesticide tax), increase in the carrying capacity characterising the wild bee populations, or a 
change in the type of insecticides used on the farm. Perhaps surprising, a decrease in the price of 
commercial bees can under some conditions help the wild bee population indirectly by allowing 
farmers to reduce their reliance on pesticides. 


II. MODELLING FRAMEWORK 

The model describes a single farm and its surrounding ecosystem supporting wild bees. Farmer 
returns consist of two components, output (yield) and costs. Output is assumed to follow a Cobb- 
Douglas production function with two inputs: pollination services and pest control. Wild bees 
together with commercial bees (if present) provide pollination services for the agricultural pro- 















4 


duction on the farm, whereas pest eontrol is aehieved by applieation of pestieides. 

A. Output 

We identify the following key eeonomie and eeologieal faetors that need to be ineluded in the 
model: (i) Both pollination and pest eontrol are essential inputs; the output is zero if either of the 
inputs is zero, (ii) Inereasing pestieide and pollination inputs inereases produetion, but only up to 
a eertain limit, (iii) Wild and eommereial bees are substitutes, but they ean perform differently for 
different erops. (iv) Pestieides affeet the population of both wild and eommereial bees, although 
the farmer ean balanee the effeet on eommereial bees by inereasing their supply, (v) Commereial 
bees, when they are present at the farm, are assumed not to eompete direetly with wild bees in 
terms of foraging for food.^ 

Thus, we assume a general form for the produetion term 

X 2 ) = A (B(y, z, Xi)f (P(xi)f , (1) 

where B is a funetion representing pollination serviees and P is a funetion representing pest 
eontrol. The level of pestieide applieation is denoted by xi, the potential for pollination by wild 
bees is denoted by y, and pollination by eommereial bees is denoted by 2 ;. a is the output elastieity 
of pollination serviees and /? is the output elastieity of pest eontrol. Further, A represents total 
faetor produetivity, and without loss of generality we assume A = 1. Funetions B(y, z, xi) and 
P(xi) are assumed to be zero if inputs are zero, P(0,0, xi) = 0 (for any value of xi) and P(0) = 0. 

The pollination serviees are assumed to inerease if the aetual population of wild bees, y, or 
aetual population of eommereial bees, 2 ;, inerease. The seeond derivatives ean be either zero 
(linear dependenee, as in the standard Cobb-Douglas model) or negative (representing deereasing 
eapaeity of bees to pollinate); in the latter ease we will assume that B reaehes an asymptote for 
large values of y and ; 2 . In order to model the effeet of pestieides on pollination, we assume that 
B{y, z, xi) is a deereasing funetion of xi. As there is not enough experimental evidenee to suggest 
a preeise dependenee, we assume a linear deerease for both the wild bees and eommereial bees, so 
that 


There is mixed evidence on this relationship in the scientific literature (Garibaldi et al. 


201 T. 
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y = K — QiXi if K — QiXi > 0 and 0 otherwise 


( 2 ) 


z = X 2 — g 2 Xi if X 2 — g 2 Xi > 0 and 0 otherwise, 
where k is the earrying eapaeity of the local wild bee population and gi, g 2 represent effeet of 
pestieides on bee population. X 2 represents the number of eommereial bees as introdueed by the 
farmer and needs to be distinguished from 2 : whieh is the aetual size of the population taking into 
aeeount death due to pestieides. Thus, the farmer buys X 2 bees, but only 2 survive and eontribute 
to pollination. X 2 together with xi are two eontrol variables in the model. 

As wild and eommereial bees are substitutes, we further simplify the funetion B{y,z,Xi) to 
depend on a linear combination of y and z,y + uz rather than on these two variables separately, u 
represents a relative efficiency of commercial bees compared to wild bees. 

Finally, we assume that P{xi) is an increasing function of xi (positive first derivative) but with 
either a zero second derivative (standard Cobb-Douglas model) or reaching an asymptote for large 
values of xi. 


B. Cost 

Wild pollination services come free (unless costly actions are taken to improve habitat), 
whereas other costs are assumed to be linear in output. The cost of pesticides is Wi per unit 
and the cost of commercial bees is W 2 per unit. We assume that the farmer minimises the cost 
function: 


c{xi,X 2 ) = WiXi + W 2 X 2 . (3) 

The minimisation is subject to satisfying a total output constraint 

q{xi,X2)=q, (4) 

where q is the exogenous target output. Given q, the levels of pesticides, xi, and commercial 
bees, X 2 , are chosen by the farmer based upon cost minimisation. 
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III. RESULTS 
A. Simplified modei 

We initially consider a simplified model in which functions B and P are assumed to be linear 
(i.e the standard Cobb-Douglas formulation). We also assume constant returns to scale and hence 
/3 = 1 — a. Finally, we neglect the effect of pesticides on commercial bees, hence g 2 = 0, and 
assume that wild and commercial bees are perfect substitutes (so u = 1). We will relax these 
simplifying assumptions later in the paper. Thus, 

q{xx,X2) = {k- giXiP X2 Tx\~°'= q, (5) 

whereas fox k — giXi < 0, 

q{xi,X 2 ) = x‘^x\~°'= q. (6) 

The farmer’s decision problem can therefore be written as 

min c{xi,X 2 ), (7) 

a;i>0,a;2>0, y=K,—giXi>{) 

q{xi,X2)=q 


1. Optimisation 

We assume that the farmer wants to achieve a target level of output, q, and hence the constraint 
becomes q (xi, X 2 ) = q (this situation can arise if farmer is contracted to achieve a certain level of 
production). A pair (xi, 2 : 2 ) represents a management strategy that she can choose to achieve this. 
The minimisation process can be represented on the (xi, 2 : 2 ) plane, with the isoquants q (xi, X 2 ) = 
q representing the output constraint. Fig. In the simple case when there is no interaction between 
pesticides and commercial bees, Eqn. 0 . the cost function is represented by a straight line with 
a negative slope, —Wi/w 2 and the intercept c/w 2 , see Fig. [^with c representing the cost. The 
procedure for optimizing, Eqn. (|7]), can then be interpreted as finding a minimal value of c such 
that the straight line still crosses the isoquant line corresponding to the given value of q, see Eig. 

The optimum value of c corresponds to the straight line that is tangent to the isoquant line; if 
more than one such line exists (representing local minima), the one with the smallest value of C is 
chosen. Eig. [^gives three examples of such lines for different values of q. 
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As the cost is linear in xi and X 2 , following standard microeconomic theory (Gravelle & Rees 


20041, this procedure can be recast in terms of the following two conditions that can be used to find 


the conditional factor demands xi and f 2 that minimise the costs of producing q units of output: 


dqixi,X2) 

dxi _ 

dq(xi,X2) 

8x2 

q{xi,X2) = q. 

The former optimality condition simply states that, in the optimum, the marginal rate of tech¬ 
nical substitution is exactly equal to the ratio of the factor prices of using xi and X 2 ', the latter is 
just the output constraint. Note that although the output function is continuous, due to the interac¬ 
tion between the two inputs it is non-monotonic and defined differently in different regions of the 
{xi,X 2 ) plane - see Fig. [^. In particular, three regions can be distinguished. Firstly, xi < k/qi 
and wild bees are present; secondly, xi > k/qi and wild bees are locally extinct. Finally, in 
a special case when X 2 = 0, pollination relies exclusively on wild bees. These three regions 
are characterised by different optimisation criteria and we first discuss each one separately. We 
subsequently compare the results to identify the strategy associated with the minimum of c. 


(8a) 

(8b) 


2. Regions 


a. Wild and commercial bees 
obtain 


For xi > 0 and X 2 > 0 the equations (8a) and (8b) can be solved to 


xi = 


where u = {1 — a)/a and 


A a "^1 + ^25' - 

X 2 = Aj*—- r - q-K 

(1 - a)w2 

y = K-giA^q. 


(9) 


Ai = 


l^W2 


( 10 ) 


wi + W2g 

The optimal values of both the pesticide use, xi, and the commercial bees use, X 2 , are both linear 
functions of the target output, q. Note that 0:2 > 0 if 


7 > 7c 


1 + 71^1 ’ 


( 11 ) 
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Commercial 


Pesticide 


FIG. 1 Farm output as a function of xi and X 2 , with values on contour lines giving g(xi, X 2 ). Straight lines 
correspond to wixi + W 2 X 2 = const. Isoquant lines and the corresponding minimum cost lines are shown 
for q = 0.8 (black), q = 1.2 (red) and q = 1.6 (blue). Vertical line corresponds to y = k — gxi = 0 and 
wild bees are extinct to the right of it. On the right, output is shown in a perspective plot. Other parameters: 
K = 1, <7 = 1, tci = 1, tC 2 = 1 and a = 1/2. 


h. Commercial bees only By assuming that y = k — giXi = 0 we can eliminate X 2 from the output 
equation and obtain the cost equation in terms of xi only, 


c{xi) = WiXi + W 2 ct^^X 


1 


( 12 ) 


Denote 



(13) 


Then, 


xi = A^q 


X2 = A^-^q 


( 14 ) 


y = 0. 

Again, the optimal values Xi and X 2 are linear functions of q? 


^ Note that, unlike in the previous section, i 2 > 0 regardless of q. 
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c. No commercial bees If X 2 = 0 then the output equation becomes 

q{xi,0) = {k - giXi)'^ 


(15) 


which can be solved to find the corresponding value of Xi. However, we can only perform an 
optimisation if this equation has more than one solution, otherwise the value of xi is completely 
determined by the output level, q. Unfortunately, no analytical solution can be found in a general 
case, but for a = 1/2 the equation is quadratic and has two solutions if q < q-m, one solution if 


q = qm, and no solutions if g > If equation (151 has two solutions, they are 


^ K - ^g\q^ 

xf = — ± -, 

' 2gi 2g 

both of which are positive. As the cost in this case is c(a;i, 0) 
solutions, Xi is optimal, hence 


(16) 

wiXi, the smaller one of these 


K - 4gig2 

“ 2gi 2g, 

£2 = 0 (17) 


y = 


K - 4gig^ 

2 ^ 2 


In this section we explore how the optimal management options change as the target output g 
increases. We show that there are three critical levels of g at which the behaviour changes: qc is 
a level at which commercial bees become economically viable, qm > qc is the maximum output 
achievable without commercial bees, and qe> qm^ qc is the output level at which the optimal use 
of pesticides leads to local extinction of wild bees. We subsequently discuss how these threshold 
levels depend on pesticide and commercial bee prices, wi and W 2 , and on the carrying capacity, k. 
We also discuss potential strategics that a social planner can use to shift the system from a state 
in which wild bees are locally extinct to the state in which they can survive. Finally, we discuss 
extensions to the model. 


B. Comparative statics 

In the previous section we have shown that the optimal strategy {xi,X 2 ) is different under 
different assumptions about the values of y and X 2 . In particular, we have identified three regions: 
(i) Region 1: X 2 = 0, (ii) Region 2: X 2 > 0 and y > 0, and (iii) Region 3: X 2 > 0 and y = 0 - 
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see Fig. As we, a priori, do not know in which of the three regions the optimal solution would 
lie for a given value of the target output, q, we first calculate the cost at the optimum, c{xi, X 2 ), 


using formulas for all three regions. As £2 can be eliminated using equation (8b), the value of 
xi that corresponds to the lowest cost is selected as an optimal value; the corresponding pair of 
the conditional factor demands, (fi, £ 2 ), describes the optimal management strategy. Let us now 
systematically discuss the various interactions in the three different regions. 


1. Region 1: Wild bees only 


We first note that for X 2 = 0 the optimisation problem has no real solutions if q > = 

^yK‘^/{2gl) (we assume here a = 1/2). Thus, qm has an interpretation of a maximum target 
output that can be achieved by using wild bees only. Note that this threshold value depends on 
ecological parameters only and not on any economic factors. We can now state: 

Proposition 1 If pollination is provided by wild bees only, there is a maximum output that can 
be achieved. This output level is determined by ecology of wild bees and their interaction with 
pesticides. 


The mechanism for this behaviour is related to the balance between pesticide use and pollinator 
population. If the farmer wants to increase the level of output, she needs to increase the level of 
pesticide use, which in turn affects the wild bee population. For small values of q, and therefore 
for small values of £ 1 , this effect is small and so the output can be increased. However, for 
large values of £ 1 , the pollinator population is reduced to such an extent that output starts to 
decline. Eventually, when £i = k/gi, the wild bee population becomes locally extinct which 
makes agricultural production impossible (as we assume that pollination is an essential input and 
it is performed here by wild bees only). 

For low values of q, Xi = x£ [cf. equation ( [T^ ] and 0:2 = 0 is the optimal choice (see Fig. 
[^). In this case, an increase in the target output, q, is possible by increasing the pesticide use 
(see Fig. |^), which results in a gradual decrease in the wild bees population (see Fig. and 
the resulting pollination services. Fig. |^). Equation (18) below shows the long run marginal cost 
(EMC), which increases non-linearly with the target output, q (see Eig. |^. That is, the farmer will 
find it increasingly more costly to increase the output by one unit as q approaches qc < 


Note that the LMC becomes infinite at q = q^n- 
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FIG. 2 Cost, c(xi,X 2 ) as a function of the pesticide use, xi, with X 2 eliminated through the constraint 
equation, q{xi,X 2 ) = q- Graphs correspond to different values of the target output, q=0A (a), 0.8 (b), 1.2 
(c) and 1.6 (d) (see also Fig. [T]l. Thin (black) line represents the case with wild bees and the thick (red) 
line represents the case without wild bees (lines are extended beyond the validity intervals to illustrate the 
behaviour; the extensions are marked as broken lines). Vertical line corresponds to y = k — gxi = 0 and 
wild bees are extinct to the right of it. Cross represents the location of the optimal solution; solid and empty 
circles represent solutions of Eqn. ( [T5] ), i.e. for X 2 = 0. Other parameters: k = 1, gi = 1, wi = 1, W 2 = 1 
and a = 1/2. 
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Qc Qm Qe Qc Qm Qe 




FIG. 3 Optimum values of (a) pesticide use, xi, (b) commercial bees population, X 2 , (c) wild bees pop¬ 
ulation, y, and (d) total pollination services, y X 2 , as functions of the target output, q. Vertical lines 
correspond to the threshold values of Qc = 0.4714045, qm = 0.5 and q^ = 1.207107, respectively. Other 
parameters: k = 1, gi = 1, wi = I, W 2 = I and a = 1/2. 


^ ^ 26igq 

dq — ^g(f 

The assumption that 0:2 = 0 is valid as long as g < q^ When q > q^ then X 2 > 0, see Fig. 
03, SO there is a sharp transition at g = g^. As qc < qm, the transition oceurs before the maximum 
possible output level aehievable by using wild bees only is reaehed (see Fig. [^. Thus, g^ is not a 
good guideline for a predietion of ehanges in the behaviour of the eombined bioeeonomie system. 
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The results can be summarised as the following proposition: 

Proposition 2 As the target output level, q, approaches Qc, it is increasingly more difficult to 
increase the output by relying on wild bees only. Thus, commercial bees become an economically 
attractive option, even though the wild bees still provide sufficient pollination levels. When q = 
fc < dm is reached, introduction of commercial bees become economically optimal. 


Qc qm Qe 




Target output 


FIG. 4 The average cost of producing a unit output (a) and the LMC, dCfdq (b), as functions of the 
target output, q. Vertical lines correspond to the threshold values of fc = 0.4714045, qm = 0.5 and 
de = 1.207107, respectively. Other parameters: n = 1, gi = 1, wi = I, W 2 = I and a = 1/2. 
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2. Region 2: Wild and commercial bees 

When q > q^, wild and commercial bees coexist, see Fig. and Fig. The pesticide use 
and the commercial bee usage increase linearly with q, while the wild bee population decreases 
linearly - see Fig. and[^, respectively. The LMC in this case does not depend on q (see Fig. 
0: 

(ir HU 4- nun 

(19) 


^ ^ + ^29 

dq ^ (1 — 0 !)w 2 


A desired increase in the target output, q, is achieved by an increase in the use of pesticides and 
the corresponding increase in the total pollination services (see Fig. [^). The increase is made 
possible by the usage of commercial bees; the total pollination levels increase. Fig. |^, although 
the wild bee population decreases (see Fig. |^). 


3. Region 3: Commercial bees only 

For values of q corresponding to Region 2, the optimal costs calculated by equations Q and 
( [T^ are similar, but c(a;i) reaches lower values if y > 0 (Region 2) - see Fig. |^. As q increases, 
both curves shift upwards, but at different rates (c/. Fig. with Figs. ^ and[^). At q = q^ 
the minimum costs using equations Q and ( [T4] ) are the same. Thus, for q = qe the solution of 
the optimality problem is not unique; selection of two different combinations of {xi, £ 2 ) (both of 
which satisfy q{xi,X 2 ) = q) results in identical cost values. This transition can be also seen in 
Fig. where the straight line corresponding to the cost function touches the isoquant line at two 
places. The threshold value can be found by equating the optimal costs calculated from Q and 
0 : 


= 


^1^2(1 — tt) 


{wi + W2g)A'^ - (1 - (^)w2 (wiAf + W2A2 “) 


( 20 ) 


For q > the optimum given by (14) is lower. This solution corresponds to Xi > n/gi and 


so to j/ = 0, resulting in the local extinction of wild bees. To understand this transition, we need 
to look at output levels which depend on both pesticide use and pollination services. At q = q^ 
the management strategy undergoes a substantial shift. Instead of relying largely on the increase 
in the pollination services (see Fig. |^), the farmer switches to a high use of pesticides (see Fig. 
[^) and a lower use of commercial bees (see Fig. [^). Note that the further increase in the target 
output will be achieved by increasing the pesticide use more than by increasing the commercial 
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bee population; compare the slopes in Figs. [3^ andjsj) for q < Qe and for q > q^.^ 

The shift in the management strategy causes the wild bees population to collapse (see Fig. |^). 
The total pollination services go down at g = q^, but continue to increase afterwards. Again, the 
LMC does not depend on q and for the parameters used here is lower than in Region 2: 

f)c 

- = w^A^ + W2Al--. (21) 

This adds another incentive for farmers to switch to the new management strategy - leading to 
extinction of wild bees - as this will not only allow to increase their output, but also to lower its 
marginal costs. 

Proposition 3 A high target output and a lower LMC can be achieved by increasing the reliance 
on pesticide use rather than pollination, leading to the local extinction of wild bees. 

C. Extending the model 

In this section we return to the general model and discuss how relaxing the simplifying assump¬ 
tions affects the results. In particular, we consider three extensions: the differential pollination po¬ 
tential of wild and commercial bees (assumption (iii)), the influence of pesticides on commercial 
bee population (assumption (iv)), and the asymptotic behaviour of the production function with 
respect to xi and X 2 (assumption (ii)). Relaxation of assumptions (i) and (v) exceeds the scope of 
this paper. 

1. Differential pollination 

First, we analyse the effect of m 7 ^ 1, i.e. different ability of wild and commercial bees to 
pollinate the crop, on the behaviour of the simplified model. In this case, the total pollination 
potential depends on ux 2 instead of X 2 . Defining X 2 = UX 2 and considering X 2 as the control 
variable leads to analogous results as for the simplified model, except that the cost function is now 

WiXi + W 2 X 2 = WiXi H- X 2 , (22) 

u 

This can also be shown by noting that Ai < A 2 , hence the slope of xi as a function of g is larger in Region 3 than 
in Region 2. 
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i.e. corresponds to a rescaled price for commercial bees, W 2 . All the results from above apply 
with the appropriate scaling of W 2 . Thus, 

Proposition 4 If commercial bees are more efficient in pollinating the particular crop than wild 
bees, M S> 1, they are more likely to be introduced (lower qj, if the unit price is the same. 

2. Effect of pesticides on commercial bees 

In the analysis above we assumed that 5 f 2 = 0 so the ability of commercial bees to pollinate 
crop is not affected by pesticides. If 5(2 7 ^ 0 then z = X 2 — 92 X 1 if X 2 — 92 X 1 > 0 and 0 otherwise. 
We will assume that the effect on commercial bees is smaller than on wild bees so that g 2 ^ di- 
Equations Q and Q can be simplified by expressing the costs in terms of z rather than X 2 , 

wixi + W 2 X 2 = wixi + W2{z + g2Xi) = (wi + g2W2)xi + W 2 Z . ( 23 ) 

Assume first that k — giXi > 0 and X 2 — g 2 Xi > 0, i.e. the pesticide use, xi, is small. Then, 

q{xi,z) = (k- giXi +z)°‘xl~°'= q. ( 24 ) 

This is the same problem as our simplified model, but with 2 ; being the control variable instead 
of X 2 , and with the price for pesticides modified by the addition if the g 2 'W 2 term. This has a 
simple interpretation: the farmer aims to control the effective population of commercial bees by 
appropriately increasing the purchased stock X 2 , with the increase given by g 2 Xi. This results in 
an increased effective price for pesticides, as their effect on commercial bees needs to be taken 
into account. 

If the pesticide use is increased so that n—giXi < 0 but X 2 —g 2 Xi > 0 (assuming here g 2 gi), 
then 


q{xi,z) = z'^xl “ = g. 


(25) 


This is again the same problem as for our simplified model. Proceeding as above, we obtain 


A2 


UW2 

Wi + W2g2 ' 


(26) 
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Then, 

Xi = A^q 
z = q 

(27) 

X2 = (A2"“ + ^2Af) g = 

^ = 0 , 

which proves that ;§ will always be positive (i.e. the use of pesticides and the input of commer¬ 
cial bees will always be adjusted so that pollination is supplied). The total cost, however, 

(wi Ag + te 2 Ag"" + w ; 25'2 A 2 ) q • (28) 

is now increased by the cost of purchasing additional commercial bees, W 2 g 2 A 2 , as compared 
to the case when ^2 = 0. Thus, 

Proposition 5 If commercial bees are affected by pesticides, behaviour is qualitatively the same 
as above, but there is an associated cost of purchasing additional bees to compensate this effect. 
This can be described by an increase in an effective unit price of pesticides. 

3. Asymptotic behaviour 

In the simplified model considered above, equations Q-Q, farm output is an increasing func¬ 
tion of the pollinator population size (both wild and commercial), but without an asymptote (see 
Fig. §. Thus, by increasing the commercial bee density one can make output to be arbitrarily 
large. In reality, the functional relationship between pollinator population size and the delivery of 
pollination services is likely to be an asymptotic type of function, so that output cannot increase 
unlimitedly. Similar features characterises pesticide use. Adding nonlinear functional forms for y, 
xi and X 2 as in equation Q significantly complicates the model and, therefore, we only consider 
Region 2 which corresponds to a high level of pesticide use and a high commercial bees density 
(and to ?/ = 0). This is the range where the saturation effect is most likely to occur; this assumption 
will be relaxed in future work. For simplicity we also assume a = /3 = 1/2. 

Given this setting, let us consider the minimization of the cost c(a;i, X 2 ) = W 1 X 1 +W 2 X 2 subject 
to the following constraint: 
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Commercial bees 
Pesticides 


FIG. 5 Functional forms for the output dependence on pollination (here on commercial bees only), for the 
standard Cobb-Douglas function, Eqn. Q with rj = 0 (thin line), and for the nonlinear form with t] = 0.1 
(dotted line), and with rj = 0.5 (thick line). The graphs are the same for pesticide use {p = 0, p = 0.1, 
p = 0.5, respectively). Horizontal lines show the respective asymptotes when applicable, also indicated in 
the margin, a = /3 = 1/2. 


q{xi,X2 ) = {B{x2 ))^^‘^ {P{xi)f^‘^ = g, (29) 

where it is assumed that y = 0. Firstly, assume that the pestieide funetional form is linear, 
P{xi) = xi (the Cobb-Douglas form still applies so that seeond derivative of q with respeet 
to a: 1 is negative), but the pollination produetion funetion tends asymptotieally to a limit aeeording 
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to: 


B{x2 ) = 


X 2 


(30) 


1 + riX2 ’ 

where l/r/ is the asymptotie value for the pollination serviee (see Fig. [^. The optimal values read 


as: 


Xl 


£2 



(31) 


c(xi,X2) = 2^wiW2g + rjwig^ . 

Note that given a = 1/2, v beeomes equal to 1 and A 2 is simplified to W 2 /W 1 {cf. Eqn. 
( [T3| )) so the above result is eompatible with equation ( [T4| ) if 77 = 0 . Thus, £2 is unaffeeted by the 
nonlinearity in the pollination produetion funetion. However, the pestieide use, £i must inerease to 
offset the relative ineffieieney of pollination. Consequently, the LMC now inereases as a funetion 
of the target output g: 


Oc 

— = 2y/wiW2 + 2T]Wig. 


(32) 


Seeondly, assume that the pollination produetion funetion is not limited, so B{x 2 ) = X 2 (the 
Cobb-Douglas form still applies so that second derivative of g with respect to X 2 is negative), but 
that the pesticide efficiency is an asymptotic function of Xi: 


P(xi) = 


Xl 


(33) 


1 + pxi ’ 

where 1/p is the asymptotic value for the pesticide efficiency (see Fig. |^. The optimal values are 


then: 


£1 

£2 



(34) 


c(fi,X2) = 2^/w/w^g + pW2q^ . 

As before, equation (34) is compatible with equation ( [T4] ) if p = 0. The solution mirrors the previ¬ 
ous case in that the optimal use of pesticides is unaffected but that the commercial bee population 
must increase to offset the inefficiency in pest control. The cost, again, increases and the additional 
term is a quadratic function of g; the LMC increases with g\ 


dc 

dg 


2y/wiW2 + 2pw2g. 


(35) 
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Finally, if both pollination and pesticide production functions have asymptotes, i.e. rj > 0 and 
p > 0, then: 



q 

c{xi,X 2 ) = - -— {2y/WiW2 + (wiP + W 2 p)q) . 

1 — rjpq^ 

In this case it is not possible to satisfy an arbitrarily increasing target output and the maximum 
possible q is given by 1/ This is due to the fact that one cannot offset a saturation in one 
factor (say, Xi) by increasing the other (say, X 2 ). Both xi and X 2 need to be set very high, resulting 
in a disproportional cost increase.^ This leads us to contend: 

Proposition 6 Introduction of saturating functional forms for xi or X 2 does not change the results 
qualitatively for Region 2, but increases the cost and makes LMC a nonlinear function of q. If 
the contribution of both pesticide use and pollination services is of an asymptotic form, there is a 
maximum output that can be achieved. 

We cannot obtain the thresholds qc and (f. analytically in this case, but one can draw some 
general conclusions by observing that higher levels of xi (and X 2 ) are needed to produce the same 
output level, if p > 0 (p > 0). This will cause wild bees to become extinct at lower levels of 
the target output, q, as compared to the standard Cobb-Douglas equation. We therefore expect the 
threshold values qc and qe to be lower in this case. 

IV. DISCUSSION 

In the paper we studied the dependence of cost-minimising management strategies on the target 
farm output, q. Thus, given the output target (for example set in the contract specifying delivery 
to a supermarket chain), the farmer will chose a certain strategy, minimising the private costs of 
output. We have shown that depending on the q, different strategies emerge, in which it is more 
profitable for the farmer to either refrain from using commercial bees when q < qc, or to use them 
if g > dc- Use of commercial bees allows the farmer to move beyond yield levels that can be 


^ We do not explicitly state the LMC in this case, but it tends to infinity as g —> 1 /y/rfp. 
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achieved by natural pollination (g > g^)- This is aehieved by inereasing both pestieide use (Fig. 
[^) and the total pollination levels (Fig. |^). 

Interestingly, we also show that there are two eompeting strategies that the farmers ean use 
to aehieve an output exeeeding q^. a low-pestieide, high-pollination strategy or a high-pestieide, 
low-pollination strategy (c/. Fig. and Fig. |^). There are good reasons to move to the latter 
management strategy, as it leads to a lower long run average eost (LAC) of produeing a unit output 
(see Fig. |^) and a lower long run marginal eost (LMC) if g > ge (see Fig. |^). However, this 
strategy also leads to loeal extinetion of wild bees if g > ge holds. 

Loeal extinetion oeeurs beeause the link between the wild bee population and output is broken, 
and the farmer does not get a signal that wild bees are deelining. This is eaused by the availability 
of eommereial bees. However, we show that the eeologieal effeet of introduetion of eommereial 
bees depends on the intensity of produetion. In the low-intensity produetion situation, the intro¬ 
duetion of eommereial bees leads to an inerease in the use of pestieides and therefore to a deerease 
in the wild bees populations (see Fig. [^. 

We also note that the transition at g = ge is an abrupt one. The wild bees population might 
be relatively low but healthy (ea. 20% of the earrying eapaeity for parameters used in this paper 
- Fig. 1^) for g smaller but elose to ge. However, for g > ge, farmers are likely to switeh to a 
high-pestieide strategy and the population of wild bees will then likely beeome extinet. Thus, a 
small ehange in the target output, eaused for example by a surge in soft fruit priees, ean make the 
high-pestieide use strategy eeonomieally attraetive, leading to a dramatie deeline in the wild bees 
population. 

The threshold output ge is therefore very important from a poliey perspeetive. The soeial plan¬ 
ner faeed with a system in whieh wild bees are loeally extinet, might want to formulate polieies 
leading to their re-establishment. The simplest poliey is to eneourage farmers to lower their target 
output, g below the threshold value, q^- Alternatively, the soeial planner might want to inerease 
the threshold value, q^. This ean be aehieved by: inereasing the priee of pestieides (inereasing 
wif, eneouraging farmers to stimulate wild bee population (inereasing k), switehing to alterna¬ 
tive pestieides (deereasing g), and/or deereasing the priee of eommereial bees (deereasing W 2 ). 


Although the formula for ge is known [see equation (20)], and it depends on both eeologieal {k, g) 
and eeonomie {wi, W 2 , a) faetors, it is not elear whether it ean reliably be estimated in praetiee. 


' This could be done by imposing a tax on pesticide use, for instance. 
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This means that the threshold might be diffieult to prediet in advanee and the transition between 


management strategies might be diffieult to ineentivise (Taylor 2009). 

While we have foeused on pestieide use, the switeh ean also be eaused by other faetors. For 
example, if the earrying eapaeity of wild bees, n, deereases due to a reduetion in the size or quality 
of wild bee-friendly habitat, the threshold values Qc and qe will deerease (see equations ( [TT] ) and 
([20|)). If the farmer still wants to attain the same target output, she will need to ehange the strategy, 
depending on the eombination of the target output and the pestieide and eommereial bees priees. 
If the target output is low, under eurrent eonditions {n = 1), the farmer does not need to use 
eommereial bees (see strategy X in Fig. |^. 

A deerease in the earrying eapaeity, for example triggered by bad weather eonditions or ehanges 
in land use eausing a reduetion in foraging or nesting areas for wild bees, leads to a shift in the 
optimal farm strategy (see Fig. - point X now lies in Region 2 demareated by the broken lines. 
Thus, the farmer might feel that the ehange in environmental eonditions forees her to introduee 
eommereial bees. 

If the level of output is high and the farmer already uses eommereial bees - strategy Y in Fig. 
1^- the reduetion in the earrying eapaeity means switehing to a high-pestieide use (as point Y now 
lies in Region 3) with an assoeiated loeal extinetion in the wild bee population. 

As noted above, the ehange in the wild bee population as the system moves from Region 2 
to Region 3 is very abrupt. For instanee, we might be dealing with what looks like a farm with 
a healthy population of wild bees (although smaller than their earrying eapaeity) for one set of 
environmental eonditions k. A small ehange in the environmental eonditions results in a limited 
loss in their ability to pollinate whieh in turns leads to the farmer modifying her management prae- 
tiees to keep up with demand for agrieultural produee. This modifieation results in a rapid loeal 
extinetion of wild bees as the threshold, qe is erossed. The importanee of maintenanee of a eon- 
dueive environment for wild bees has also been emphasised by |Keitt| ( |2009| ), who eoneluded that 
bee populations would deeline abruptly if environmental eonditions (in this ease, habitat) deteri¬ 
orated below a eertain threshold. The results of this paper support this eonelusion, as inereasing 
the earrying eapaeity of the surrounding area would allow higher outputs without exeeeding the 
threshold leading to population extinetion. 

The partieular funetional forms are ehosen in this paper for their general applieability to a wide 
range of agri-eeologieal problems as well as for their simplieity; this applies to both the linear 
dependenee of B on xi and to the Cobb-Douglas produetion funetion. The latter has been ehosen 
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1^2 = 0.5 



FIG. 6 The effect of changing k on the ecological outcome under the optimal management strategy as 
a function of the marginal pesticide cost, wi and the target output, g; k = 1 (solid lines) is reduced to 
K = 0.5 (broken lines). Lines represent the threshold budget values, ijc (black line) and qe (red line). Other 
parameters: gi = I and a = 1/2; m 2 = 0.5. X and Y strategies are discussed in the text. 


over alternative forms to model pollination serviees beeause it refleets a eomplete relianee of the 
produetion system on wild or eommereial bees. This assumption will be aeeurate for many erops 
with a high dependenee on pollinators, sueh as many berries and orehard fruits ([Hanley et al. 


20151. In partieular, if y and 2 ; are both zero, the output is zero, independently of Xi. Although 


in our paper xi represents the use of pestieides, it ean also be interpreted as any agrieultural 
praetiee that (i) is essential for the output generation, and (ii) affeets wild bee populations. This 
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generalisation strengthens the ease for the partieular funetional form. Finally, we assumed that 
eommereial bees do not direetly eompete with wild bees. This assumption ean be relaxed and we 
expeet this to make the population of wild bees more fragile. This extension would require a more 
realistie population model and exeeeds the seope of this paper. 

In our model the ehanges are fully reversible (the wild bees population reaets immediately to 
ehanges in xi), but in reality sueh shifts are likely to be irreversible. This will be problematie 
if the strategy taken then beeomes uneeonomie, for example due to an inerease in the priee of 
pestieides, or eommereial bees. Sueh a switeh has already oeeurred in the apple growing region of 
Siehuan, China, where human pollinators were used as substitutes, allowing a high pestieide, low 
habitat strategy to eontinue. When human pollination beeame too expensive, the only option for 
farmers was to leave the market altogether, and diseontinue apple produetion. When deelines in 
wild eapital sueh as wild pollinators are irreversible, and there is uneertainty over whieh souree of 
eapital will be most benefieial in the future, there is a value to maintaining the natural eapital for 
future use ( Arrow & Fisher[ 1974 Kassar & Lasserre 20041. This “option” value is an ineentive 
for eonserving wild pollinators, and will be positive even if there are no immediate advantages of 
supporting wild pollinators. 

The wild bee population modelled here will often be made up of multiple populations of bee 
and non-bee pollinators (sueh as hover flies). The presenee of multiple pollinator groups eould 
buffer the system to extinetion; the relative toleranee of pollinator networks to extinetion has 
been shown by [Kaiser-Bunbury et al (2010); Memmott et al. (2004). However these studies 
do not assume that threats to the different populations are eorrelated. While different pollinators 
groups may respond in slightly different ways to external pressure sueh as pestieide use, the effeets 
are likely to be negative on all groups, and may be stronger on non-bee pollinators as these are 
smaller ( Goulsonl |2013| ). The model diseussed in this paper is unique in its inelusion of a ehronie 
threat to pollinators (pestieide use), whieh is likely to affeet all pollinator groups. The benefit 
of maintaining multiple groups of eeosystem serviee providers as insuranee against a fluetuating 
environmental was diseussed by [Baumgartner (2007). The problem eonsidered here differs as we 
eonsider a threat whieh is likely to be detrimental on the whole pollinator eommunity, means that 
holding diverse pollinators will not be benefieial, however maintaining both eommereial and wild 
bees will be valuable as options for pollination provision in the future. 

In our model farmers aet myopieally. We also assume that wild bees respond instantaneously 
to the ehanges in the management. In effeet, the model demonstrates behaviour that would be 
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observed if the farmer can make planning decisions once and see the impact of those in the future 
without the chance of adaptation. This may well be realistic, as farmers are unlikely to be able to 
detect small changes in wild bee populations from year to year, but will notice dramatic decreases 
in pollination services over longer periods of time. However, the model, and in particular its agent- 
based extension, can be generalised to include different planning horizons for farmers as well as 


the long-term dynamics of bees, for example in the form of a dynamic model of Khoury et al. 
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